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Abstract 

In this paper, theoretically, the characteristics matrix method is employed to investigate and compare 
the properties of the band gaps of the one-dimensional ternary and binary lossy photonic crystals which 
are composed of double-negative and -positive materials. This study shows that by varying the angle of 
incidence, the band gaps for TM and TE waves behave differently in both ternary and binary lossy structures. 
The results demonstrate that by increasing the angle of incidence for the TE wave, the width and the depth 
of zero-n, zero-/i, and Bragg gap increase in both ternary and binary structures. On the other hand, the 
enhancement of the angle of incidence for the TM wave, contributes to reduction of the width and the depth 
of the zero-ri and Bragg gaps, and they finally disappear for incidence angles greater than 50° and 60° for 
the binary structure, and 40° and 45° for the ternary structures, respectively. In addition, the details of the 
edges of the band gaps variations as a function of incidence angle for both structures are studied. 
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I. INTRODUCTION 



Photonic crystals (PCs) are artificial dielectric or metallic structures in which the refractive 
index changes periodically. The periodic structure of the PCs contributes to photonic band gap 
(forbidden range of frequencies). Interference of the Bragg scatering is considered as a cause of this 
phenomena. PCs have attracted much interest due to their novel electromagnetic waves character- 
istics and important scientific and engineering applications and have received attentions by many 
researches in recent decades [l-8]. With possibility of producing metamaterials, PCs containing 
metamaterials, namely "Metamaterial Photonic Crystals" (MPCs) have been made. The techno- 
logical relevance of the MPCs is of great importance. These structures have attracted considerable 
attentiom for their various applications. Designing dielectric mirrors, effective waveguides, filters, 
and perfect lenses are clear cut examples. The metamaterials are utalized to design binary and 
ternary MPCs with layers of Double Negative (DNG) and Double Positive (DPS) materials. In 
these structures, various band gaps can be seen in the transmission spectrum. When the average 
refractive index of the MPC is equal to zero, the band gap is called zero-n gap jjj. The properties 
of the zero-n gap of one-dimensional (ID) binary MPC that composed of the two layers of DNG 
and DPS materials were investigated in recent years [9l4l6t|. In ID binary PMC, in addition to the 
zero-n and Bragg gaps, we also have zero-u and zero-e gaps for TE and TM waves respectively, 
which appear for the oblique incidence angles. The results have been also reported by some authors 



lfHIS]. 



The novel idea of this research is to include the loss factor which is an unavoidable consequence 
of Double Negative materials to investigate optical properties of ID binary and ternary MPCs. In 
this regard the transmission of electromagnetic waves through a ID ternary lossy PC consisting of 
layers with double-negative and -positive materials is studied. The effects of the incidence angle and 
polarization on the characteristics of the zero-n, zero-u, zero-e, and Bragg gaps are investigated. 
Furthermore, the results are compared with the similar results in the ID binary lossy PC. 

The paper is organized as follows: the MPC structures design, the permittivity and permeability 
of the DNG layer, and theoretical formulation (characteristic matrix method) are described in 
Section 2, the numerical results and discussions are presented in Section 3, and the paper is 
concluded in Section 4. 
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II. STRUCTURES DESIGN AND CHARACTERISTIC MATRIX METHOD 



The ID MPC structure under study which is located in air is constituted by alternative layers 
of DNG and DPS materials, where the DNG material is dispersive and dissipative. We consider 
the ID binary MPC with periodic structure of (AB) N and the ID ternary MPC with periodic 
structure of (ABC) N , where A, B, and C denote three different materials. N is number of the 
lattice period and di, Si, and /ij (i = A,B, C) are thickness, permittivity and permeability of the 
layers, respectively. 

The permittivity and permeability of the DNG (layer ^4) in the microwave region are complex 
and are given as [9], 
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where / and 7 are frequency and damping frequency, respectively, given in GHz. Various details of 
the real parts of the permittivity and permeability of layer A, versus frequency have been discussed 
in our previous work 110 1. 

□ 

The calculation is based on the characteristic matrix method [19|], which is most effective to an- 
alyze the transmission properties of PCs. The characteristic matrix M[d] for TE wave at incidence 
angle from vacuum to a ID PC structure is given by []jj]: 



Mi 
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-1 pi sin 7, cos7i 

where k = 2 and 3 for binary and ternary PCs, respectively, ji = (uj/c) n^di cos 8i, c is speed of 
light in vacuum, 9i is the angle of refraction inside the layer with refractive index rii and pi is given 
as: 
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The characteristic matrix for N period structure is therefore [M(d)] N . The transmission coef- 
ficient of the multilayer is given by: 

2 Po 
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(6) 



where rriij(i,j = 1,2) are the matrix elements of [M(d)] N , and 

po = n cos6 ,p s = n s cos9 s (7) 

The transmissivity of the multilayer is given by: 

T = ^\t\ 2 . (8) 
Po 

The above formulations can be applied for TM wave by simple replacements of pi, po, and p s 
as follows: 

Pi = J— cos6i (9) 

V 
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III. NUMERICAL RESULTS AND DISCUSSION 

Based on the theoretical model described on the previous section, the transmission spectrum of 
the presented lossy MPC structures was calculated. In the study of the ID binary MPC, consisting 
of DNG (layer A) and DPS (layer B) materials, equations (1) and (2) are used for the permittivity, 
ea, and the permeability, ha, of layer A. Layer B is assumed to be a vacuum layer with eb = l±B = 1 
(rig = 1). The thickness of layers A and B are chosen as o\a = 6mm and ds = 12mm, respectively 



id 



and the total number of lattice periods is set as N = 16 

The transmission spectra of TE and TM polarized waves for the binary structure at various 
angles of incidence and for 7 = 0.2 x 10 -3 GHz are shown in Figs. 1(a) and 1(b), respectively. 
For oblique incidence the zero-/i and zero-e gaps appear in the transmission spectrum for TE and 



TM waves, respectively, as reported in 



13, 



181 ] . The gaps appear at the frequencies where the sign 
of the permeability or permittivity of DNG materials changes. As it is seen in Fig.l the width of 
the gaps increases as the incidence angle increases while the left edge of the gaps remain nearly 
unchanged. In addition, for TE waves the width and the dep th of the zero-n and the Bragg gap 
increase as the incidence angle increases, as reported in [9- 14] . On the contrary, for TM waves the 
width and the depth decrease as the incidence angle increases (Fig. 1(b)). Moreover, the figures 
show that the central frequency of the zero-n, as reported in our previous work and the Bragg 
gap shift to the higher frequencies as the angle of incidence increases. 

In the next part, the band gaps of the ID ternary MPC are investigated. The binary structure 
used before, is modified by introducing a third layer of Si02 (layer C) with refractive index of 



nc = 1-46 and thickness of dc = 6mm in each lattice period. The other parameters are kept 
the same as in the binary structure. The transmission spectra of TE and TM polarized waves 
for the ternary structure at various angles of incidence and for 7 = 0.2 x 10~ 3 GHz are shown in 
Figs. 2(a) and 2(b), respectively. As it is seen, the width of the zero-e and zero-// gaps increases 
by increasing the incidence angle. Although, for TE wave the width and the depth of the zero-n 
and the Bragg gaps increases as the incidence angles increases. For TM waves the behaviour is 
completely different, such that the width and the depth of the band gaps decrease when the angle 
of incidence increases. However, the central frequency of the zero-n and Bragg gaps behave the 
same as in the binary structure where they shift to the higher frequencies as the angle of incidence 
increases. 

In this part the band gaps in binary and ternary structures are compared. As it is seen from, 
Figs. 1 and 2, the width of the band gaps decreases, when another dielectric layer added to the 
lattice period of binary structure. Furthermore, the zero-n and Bragg gaps frequencies of the 
ternary structure are lower the corresponding frequencies in the binary structures; while the zero-e 
and zero-/i gaps appears in the same frequencies in both structures. 

The lower (/x,) and higher (/# ) frequencies of the zero-n gap in the binary and ternary struc- 
tures, as a function of incidence angle for TE and TM waves are shown in Figs. 3(a) and 3(b), 
respectively. As it is clearly seen, the lower frequency of the gap for the ternary structure is smaller 
than that in the binary structure for both TE and TM modes for different incidence angles. It is 
also interesting to note that the width of the gap in the ternary structure for TE wave is sensitive 
to the angle of incidence and increases as the incidence angle increases while the left edge of the 
band gap is nearly invariant. Moreover, the width of the zero-n gap of both binary and ternary 
structures is very sensitive to the incidence angle for TM wave, and the gap disappears for inci- 
dence angles greater than 50° in the binary structure as reported in [10,12]. The gap disappears 
for incidence angles greater than 40° for the ternary structure. 

The lower and the higher frequencies of the zero-// and zero-e gaps as a function of incidence 
angle, for both binary and ternary structures are shown in Figs. 4(a) and 4(b), respectively. The 
width of the zero-// and zero-e gaps increases faster for the binary structure than the ternary ones. 
As it is clearly seen, the gaps appears at the same frequency for both binary and ternary structures; 
in addition, the figures indicate that the upper edge of the gaps, so the width of the gaps, are very 
sensitive to the incidence angle while the lower edge of the gaps are nearly insensitive. 

The Jl and fjj frequencies of the Bragg gap of the binary and ternary structures, as a function 
of incidence angle for TE and TM waves are shown in Figs. 5(a) and 5(b), respectively. It is seen 
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that the Bragg gap like the zero-n gap appears in the lower frequencies for the ternary structure. 
Moreover, /x and fn shift to the higher frequencies as the incidence angle increases, and as it is 
seen from Fig. 5(b), the gap for TM polarized wave disappears for incidence angles greater than 
60° for the binary and 45° for the ternary structures. 



IV. CONCLUSION 



The numerical results show that including the loss factor in the permittivity and permeability 
of the DNG layer the changes in transmission spectrum of TE wave in ID ternary and binary 
PMCs behave similarly for different incidence angles. The zero-n and zero-//, and the Bragg gaps 
become wider and dipper as the angle of incidence increases. In addition, by enhancing the angle of 
incidence, the lower frequencies of the zero-n and zero-yLi gaps are nearly invariant in both ternary 
and binary PMCs, but the Bragg gap increases. Moreover, the upper edge of all three band gaps 
increases. Such behaviors do not clearly observed in the zero-n gap. The transmission spectrum 
for TM polarized wave is rather different. It is found that in both structures, as the angle of 
incidence increases, the width and the depth of the zero-e gap increase, while for the zero-n gap 
and for the angles of incidence greater than 40° for the ternary and 50° for the binary structures the 
zero-n gap disappears. Similar behavior is observed for the Bragg gap for the angles of incidence 
greater than 45° and 60° for binary and ternary structures, respectively. For TM waves, the lower 
frequency of the zero-n and Bragg gaps increases, but the zero-e gap decreases as the angle of 
incidence increases in both binary and ternary structures. The higher edge of the zero-e and Bragg 
gaps increases while it decreases for the zero-n gap. It is interesting to note that the difference 
between the lower and higher frequencies of all the gaps for TM wave in the binary PMC, is larger 
than that in the ternary PMC as all the figures show. To recapitulate what was said before, i- 
the width of the band gaps decreases as a dielectric layer is added to the binary structure, ii- the 
zero-n and Bragg gaps appear at the lower frequencies in the ternary structure, iii- the zero-e and 
zero-// gaps appear at the same frequencies in both structures. Finally, the results of the present 
study could be employed in designing new edge filters, waveguides, and other optical devises in 
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microwave engineering. 
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FIG. 1: Transmission spectra for the ID binary MPC structure, for different angles of incidence, with 
7 = 0.2 x 1(T 3 GHz; (a) TE (b) TM waves. 




FIG. 2: Transmission spectra for the ID ternary MPC structure, for different angles of incidence, with 
7 = 0.2 x I0" 3 GHz; (a) TE (b) TM waves. 
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FIG. 3: The lower and higher frequencies of the zero-71 gap as a function of angle of incidence for both 
binary and ternary structures; (a) TE and (b) TM waves. 
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FIG. 4: The lower and higher frequencies of the (a) zero-// and (b) zero-e gaps as a function of angle of 
incidence for both binary and ternary structures. 
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FIG. 5: The lower and higher frequencies of the Bragg gap as a function of angle of incidence for both binary 
and ternary structures; (a) TE and (b) TM waves. 
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